Helicase activity on DNA as a propagating front 
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We develop a propagating front analysis, in terms of a local probability of zipping, for the helicase 
activity of opening up a double stranded DNA (dsDNA). In a fixed- distance ensemble (conjugate 
to the fixed-force ensemble) the front separates the zipped and unzipped phases of a dsDNA and a 
drive acts locally around the front. Bounds from variational analysis and numerical estimates for 
the speed of a helicase are obtained. Different types of helicase behaviours can be distinguished by 
the nature of the drive. 



A helicase moves along a DNA and unwinds it. When- 
ever a double stranded DNA in a cell needs to be opened 
up, a helicase is summoned, be it during the semi- 
conservative replication, repair mechanism of a stalled 
process or other DNA related activities A large num- 
ber of helicases including rna-helicases have so far been 
identified from different pro- and eucaryotes. A well- 
studied bacterium like E. Coli contains at least 17 differ- 
ent helicases, though the need and the function of each 
of these are not yet clear. 

The helicase activity involves a motor action fed by 
NTP's (nucleotide triphosphate) and eventual opening 
of DNA by dissociating successive base pairs along the 
chainQ. Quantitative estimates of rates of such activ- 
ities (~ 400 base pair per second or less) for almost 
all helicases are known from in- vitro studies in solu- 
tions and more recently from single molecular experi- 
ments. Attempts to categorize these varieties of helicases 
as per their common features have led to various classifi- 
cation schemes. These are: (i) active vs passive depend- 
ing on the direct requirement of NTP for the opening; 
(ii) families and superfamilies (SF) based on the con- 
served motifs of the primary sequence; (Hi) monomeric, 
dimeric, hexameric, oligomeric depending on the number 
of units required for activity; and (iv) mode of translo- 
cation: whether it translocates on the single stands of 
DNA or on dsDNA. For example, dnaB, the main he- 
licase involved in the replication of DNA in E. Coli, is 
a hexameric, passive helicase belonging to the dnaB-likc 
family translocating on a single strand DNA0|. PcrA 
is an active, SF1, monomeric helicase translocating on 
ssDNAQ while recG is an active, SF2, monomeric, ds- 
DNA translocating helicase^. Apart from these gross 
classifications, very little information is available on the 
detailed mechanism of the helicase activity. 

Crystallographic data available for a few helicases have 
been used to model mechanisms for specific helicases. 
Though crystal structures cannot give a dynamic view, 
such proposals, attractive no doubt, are the only ones 
available so far. According to these proposals, hexam- 
eric helicase like dnaB, opens up dsDNA like a wedge 
by virtue of its motor action on ssDNA0|. A rolling 
mechanism has been proposed for dimeric helicases |3- 
In case of PcrA, the helicase activity and the motor ac- 
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FIG. 1: Schematic representation of the advocated mecha- 
nisms for (A) dnaB and (B) PcrA helicases (see text) whose 
motions are indicated by the arrows. (C) Proposed model. 
Constraints are that the DNA is open on one side but zipped 
on the other, both phases coexisting. The unzipped phase is 
preferred within the shaded box of width wh encompassing 
the Y-fork or interface or domain-wall. 



tion can be decoupled. Crystal structure, supplemented 
by biochemical evidences on mutants, favours a mecha- 
nism where the helicase moves forward on ss-DNA and 
during its sojourn a different domain of the helicase pulls 
a few bases of a strand of the ds-DNA beyond the Y-fork, 
the junction between ss and ds-DNA0]. See Fig. [I] for a 
schematic diagram. RecG is more complicated because 
it moves in opposite direction from zipped to unzipped, 
a fork reversal process forming a Holliday junction of 4 
dna strands 0. 

Our aim is to develop a generic physical picture that 
could be applicable to all the different types of heli- 
cases. Recently, a phase coexistence based mechanism 
for helicase activity has been proposed by Bhattachar- 
jee and SenoQ. A kinetic model has also been proposed 
recently^ while a random walk model was used in an 
earlier study to analyze the movement on DNA. The 
phase-coexistence mechanism is based on the unzipping 
phase transition of a ds-DNA by a force, which was first 
shown in a continuum model in Rcf. [lOj. This transi- 
tion has since then been established by exact calculations 
for lattice models in all dimensionsllll |l2l. in studies 
of dynamics [HI El, by scaling theories|l4j. in quenched 
averaged DNA |l5|. and others. The key points of the 
mechanism are the following, (a) A helicase keeps two 
single strands of a dsDNA at a separation much bigger 
than the base-pair separation of a ds-DNA. (b) The re- 
sulting fixed-distance ensemble for a dsDNA breaks it up 
into a zipped phase and an unzipped phase separated by 
a domain wall, (c) The domain wall is identified as the 
Y-fork junction. All helicases act at or near this junc- 
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FIG. 2: Plots of (A) V(p) and (B) V(cf>). Curves (a) show 
the co-existence of the zipped (p = po or <f> = 0) and the 
unzipped (p = or (f> = 1) phases. Curves (b) and (c) show 
the change near the interface or domain wall due to the drive 
(h in Eq. ||1J) that makes the zipped phase locally metastable 
or unstable. 



FIG. 3: Schematic representation of the pulse: The front is 
at location zi (filled circle at <j> = 0.5) at an arbitrary time 
t\. The effective free energy is of form (b) or (c) of Fig. [5] 
(metastable or unstable) only within z\ ± Wh- The drive or 
pulse moves with the front. The relaxation due to the elastic 
term leads to a traveling wave solution. 



tion of the ds-ss DNA . (d) the motor action of the heli- 
case leads to a shift in the position of the fixed-distance 
constraint thereby shifting the domain wall towards the 
zipped phase. Additional features are needed for effi- 
ciency, job-requirement and processivity (the distribution 
of the length unzipped before a helicase drops off). Our 
purpose in this paper is to use this coexistence hypothesis 
to develop a simple coarse grained model of the propa- 
gation of the Y-fork. 

We use the zipping probability p{z), the probability 
that the base pair at index or contour length z (measured 
along the backbone) is zipped as the basic variable. The 
unzipping transition by a force has hitherto been stud- 
ied by using polymer models. However, in the case of 
the conjugate fixed distance ensemble, the probability of 
zipping, p(z), under the imposed constraints, has been 
shown to be a useful description 0, 01- In the polymer 
approach, one introduces a local variable n(z) = 1 or 
depending on whether the base pairs at z are bound or 
not to write the pairing energy, characteristic of the DNA 
problem, as J dze(z)n(z) where e(z) is the base pairing 
energy at z. The average value (n(z)) gives the fraction 
Po of bound basepairs, and it is the parameter monitored 
in melting and unzipping transitions. This suggests the 
use of p(z) as a coarse-grained parameter to use for the 
inhomogeneous case, especially for the description of the 
interface. Since the unzipping transition is of first or- 
der nature with a co-existence on the phase boundary, 
and our interest is in the region much below the thermal 
melting point of DNA, a phenomcnological Landau-type 
hamiltonian or free energy can be used to describe the 
state of the DNA. In addition, the presence of the heli- 
case can be represented by the boundary conditions on 
two sides of the DNA (zipped on one side and unzipped 
on the other) . The effect of the motor action that pushes 
(or pulls) one phase into the other and any other desta- 
bilizing effects of the helicase are taken into account by 
additional terms involving the zipping probability p{z). 
We show the existence of a traveling wave solution (i.e. 
a propagating front of Y-fork) and then get the speed of 
the front in terms of the parameters of the hamiltonian. 



A variational principle is used to derive bounds on the 
selected velocity. 

The unzipping transition can be described by a 
Landau-type hamiltonian or free energy 



Ho = dz 



^p{ z f + V{p{z)) , (1) 



with 



v (p) = 2 r P 2 + 3 W P 3 + 4 U P 4 ' 

where K is the appropriate rigidity modulus, p(z) = 
dp{z)/dz, and the cubic term ensures a first order tran- 
sition for w < 0, u > 0. See Fig. [21 The coefficients r, w 
and u can in principle be determined from the knowledge 
of the unzipping phase boundary and po- It is often use- 
ful (and used here also) to reparametrize V{p) in terms 
of p , the width of the well at po and the barrier height 
AE. Treating Eq. Q as the meanfield free energy [13, 
the dynamics of unzipping is given by the overdamped 
equation of motion 



,-i dp(z) 
dt 



5p(z) 



(2) 



with appropriate transport coefficient T. Since we are 
away from critical points and interested in nonequilib- 
rium propagation problem, we may ignore noise terms in 
the equation of motion. Stochastic terms would also be 
required to describe processivity, which we do not con- 
sider in this paper. 

For a long chain, with the boundary conditions p(z) — > 
po (probability in the zipped phase) as z — > +00 and 
p(z) ^ as z - 00 at phase coexistence, the vari- 
ation of the equilibrium probability of zipping is de- 
scribed by the equation Kp(z) = V'(j>), with prime 
(dot) denoting derivative with respect to the argument 
(z). For the assigned boundary condition, there is a 
domain wall solution located at an arbitrarily chosen 
z = with a profile f£ dp/ y / 2V(p)/K = z, and p(z) ap- 
proaches the bulk value exponentially in the zipped phase 
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as po — p(z) ~ cxp(— V"(po) | z | /yK). The energy of 
a wall of width W, in terms of the barrier height AE 

2 

(Fig. is E = Kjfir + AE W, which on minimization 
gives W = (Kp^/AE) 1 / 2 . Measurements of W, and of 
AE (from the decay rate of bubbles in the bound state) 
would give an estimate of K . This equilibrium situation 
can be obtained by keeping the helicase static on dsDNA. 
e.g. by denying ATP in in vitro experiments. Since such 
configurations can now be prepared jl8j|. detailed charac- 
terization of the wall can be done experimentally. 

In order to incorporate the effect of the motion of the 
helicase, we introduce a moving perturbation that tends 
to destabilize the domain wall at its current location. A 
time dependent perturbation is introduced in the equa- 
tion of motion, Eq. or, cquivalently, in the Hamilto- 
nian of Eq. I|T)1. that favours the unzipped state over a 
region of width Wh around the domain wall, maintaining 
co-existence elsewhere. See Figs. |21 and |3] The crystal 
structure of PcrA0| suggests Wh ~ 20 bases. The equa- 
tion of motion, augmented by a "drive" term, is now 
given by 

= K^-rp + ™P 2 -up s -h(z,t)V{(p), 



with h{z,t) = U 



z — ct 
w h 



(3) 



This is equivalent to adding a term J dzh(z,t)Vi(p) in 
the Hamiltonian, Eq. QJ, such that the drive favours the 
unzipped region (p = 0) in a region of width Wh around 
z = ct, with the front position at t = chosen as origin 
z = 0. In Eq. V\{p) should have the right form 

for V(p) + Vi(p) to favour the unzipped phase. Since 
the helicase works only near the interface or the front, 
Li (x) is a short range function. For simplicity, we choose 
U{x) = Au for | x |< 1 and zero otherwise. The "drive" 
is attached to the front and moves to the zipped side with 
a speed c which is to be determined self-consistently so 
that the front also moves with the same speed. 

The role of the drive term is to disturb the coexistence 
between p = (unzipped) and p = pa (zipped phase). 
By a transformation of variables, like p = po — <f> and 
rescaling, we recast Eq. (J2J in a more standard and sym- 
metrical form (by choosing V±(p)) 



d<f>(z) 
dt 



-—2 +f{4>), where 
az z 



(4) 



1 



+ h{z,t) 4> + 



+ h(z,t) 3 , 



and, for brevity, same notation h is used for the drive. 
Another choice could have been /((/>) = 0(1 — 4>){<j> — 
h — h(z,t)) which is identical to Eq. |J3J upto a scale 
transformation if and only if h = constant. In any case 
no fundamental difference is expected among the various 
possible choices. Eq. allows the zipped phase with 



= and the unzipped phase with = 1 for all h. These 
two phases coexist at h = 0, while <f> = is metastable 
for < h < 1/3 and unstable for h > 1/3. See Fig. □ 
The symmetrical form helps in identifying the location 
of the front by <p = 0.5 and the drive h(z,t) is operative 
only around that point. 

Assuming that the front propagates with a velocity c, 
i.e. <j)(z, t) = <j)(z — ct), We can use the comoving frame 
with coordinate £ = z — ct to rewrite Eq. |@Jas 



d 2 cj> d(j>{z) 



+ c 



d ~ 



d ~ 

with -t-V(4>) 



-fit) 
(5) 

which has a mechanical analogy of a particle moving in 
a potential —V(4>) under friction (taking £ as a time 
like variable). This analogy immediately tells us (using 
first integral or energy conservation) that to satisfy the 
boundary conditions at £ = ±oo when h = 0, one must 
have c = 0. In other words, there is no propagating so- 
lution as it should be in the case of phase coexistence. 
A propagating solution ensues if the drive h is not zero, 
with c = c(Au,Wf l ). 

The speed of propagation c has to be insensitive to 
width Wh of the pulse if Wh » the width of the front 
or interface. In that large Wh limit, c should be the 
speed of a propagating front for a uniformly metastable 
(i.e. h < 1/3) or unstable (i.e. h > 1/3) case. We 
recollect the relevant results for the uniform situation 
with h(z,t) = constant, (i) There is a pushed to pulled 
transition[17j in the propagation of the front at h = 4/3. 
(ii) Beyond h = 4/3 the velocity is determined by the 
linearized equation of motion while the full nonlincar- 
ity is important for h < 4/3. (in) For the metastable 
case {h < 1/3) any initial condition <po(z) = <fi(z,t = 0) 
rapidly evolves to a steady shape with a velocity c' (h) = 
3/i/^/2[(2/3) + h], approaching the steady speed expo- 
nentially fast in time, (iv) For the unstable case, a sharp 
interface (say a sharp step at t = 0) also evolves to this 
"pushed" front so long as h < 4/3. However a flatter in- 
terface would maintain its flatness and move with a speed 
determined by the initial flatness, (v) In the pulled limit, 
(h > 4/3), the asymptotic speed is the linearized speed 
c*(h) = 2y/i — (1/3) which is reached algebraically in 
time provided the initial condition is sharp (e.g. a step 
function) . 

With a pulse, one may use the particle mechanics anal- 
ogy that a particle is in one of the peaks of the equal 
height double peaked potential —V{4>) at £ = — oo and 
then at finite time it gets a push (energy input in par- 
ticle mechanics) which should be sufficient to overcome 
frictional loss and reach the other peak at £ = +oo. This 
will satisfy the boundary conditions at £ = ±oo. A mov- 
ing front is therefore possible. In other words the mov- 
ing front originates from the "elastic term" that tries to 
spread out the change in (f> in the front region. 
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FIG. 4: Positions of the front or domain-wall for Am = 0.2 
(see below Eq. |5J at bulk coexistence are shown by plotting 
}>{z) vs z. The curves are at times 0, 200, 400, 600, 800, 1000 
in arbitrary time units and are obtained from a discretized 
version of Eq. J?) . 

For a quantitative analysis of the speed, we use a 
variational principle |l9j. If the equation of motion ad- 
mits, as we verify numerically below, a monotonic front 
4>{z,t) = <?(£), then the inverse mapping can be used to 
get £ from q with < q < 1. There is an inequality, 

c 2 > 2 /° f9dq , (6) 

for any positive function with —dg/dq > 0, provided 
the integrals exist. This requires /'(0) < for bistable 
/ of the form Eq. (0} . Only the metastable case is 
considered here. The unstable case with /'(0) > can 
also be treated though in a slightly different way. Taking 
g(q) = [(1 — q)/q] b and the uniform case of / i.e. h(z,t) = 
h, one gets 

c 2 > {-b 2 /3) + hb(6-b)/2. (7) 

The supremum of the lower bound at b = 9h/(2 + 3h) 
recovers the exact velocity, c\ mentioned earlier (valid 
for h < 1/3). Under the assumption of monotonicity, the 
pulse, in q-space, is at q = 1/2 and is of width Cl(u)h) 
such that fl(Wh — > oo) = 1/2. The profile q(£) gener- 
ally approaches the two limits exponentially as we have 
seen earlier and therefore Cl = (1/2) — Q ~ cxp(— Wh/w*) 
with some characteristic length w* . The pulse term con- 
tribution to the numerator of the bound of Eq. JSJ) is an 
integral of q(l— q 2 )g(q) over q e [f2, 1— ti] and the integral 
can be expressed in terms of incomplete beta functions. 
We then obtain 

-b 2 1 

c 2 > + -Aub(6-b)-A(n,b)Aub, (8) 

3 2 

where the form of A(Cl,b)(> 0) is not displayed. Taking 
the maximum of the right hand side as the best estimate 
c(Au 7 Wh) for the speed, we see that c(Au,Wh) < (the 
bulk value), as expected. Using the asymptotic behaviour 
of the incomplete beta functions in A(Cl, b), one finds that 
c(Au,w) saturates exponentially for large Wh- In the 



other limit of small Wh (equivalent to small f2 ), there 
is a linear dependence on Wh- Combining these, we may 
write c(Au, w) = c [l — aexp(— w/wo)], a form that does 
represent the numerical data very well. 

In order to determine the velocity of propagation of 
the front with a pulse, we have numerically solved a dis- 
cretized version of Eq. JSJ) for various values of the width 
u>h and magnitude of the drive Ait. A small time step is 
chosen for proper convergence but no spatial continuum 
limit has been done. We start with a sharp interface, 
4> = 1 for < z < L/2 and = for z > L/2. A se- 
quential update is done. At every time step, we allow a 
square pulse of width Wh and strength Au at the current 
location of the domain wall or front (located by 4> = 0.5). 
In all cases we observed a monotonic front. Fig. 01 shows 
the time evolution of the front for the case of a drive with 
Au = 0.2 and of zero width, Wh = (pulse at one point 
of the lattice only ). 

The variation of the speed with the width of the pulse 
is shown in Fig0 Consistent with our results from the 
variational analysis, we see that the velocity approaches 
the bulk limit for large widths and this approach is ex- 
ponential. There is a small but systematic deviation of 
the observed velocity from the exponential fit for larger 
values of Au (in the "unstable" region). Detailed anal- 
ysis of the pushed versus pulled cases will be reported 
elsewhere. 

In terms of helicases, it seems natural to associate 
dnaB type passive helicases with the metastable case 
where the motor action provides the drive that locally 
disturbs the Y-fork region. As in Ref . Q , we identify the 
domain wall or the front as the Y-fork region - the junc- 
tion of the ds-ss DNA. In the metastable case, the pushed 
dynamics has a stability against small fluctuations, the 
speed of propagation is independent of the initial con- 
ditions, and the speed approaches the steady state limit 
exponentially in time. All of these are important prop- 
erties expected of a hclicasc of type dnaB which, after 
loading on DNA, carries out the unzipping in tandem 
with the other processes during replication. 

So far as PcrA (Fig. ^) is concerned, we associate the 
overall dynamics to the unstable case. No quantitative 
experimental results are available regarding the magni- 
tude and width of the force PcrA exerts on the bases 
beyond the Y-fork. It is reasonable to assume that the 
force is meant to unzip DNA locally, and the effect of this 
pulling is to make the ds-region unstable. We then infer 
that PcrA operates in the unstable (h > 1/3) regime. 
In several mutants of PcrA (replacing a few residues by 
alanines) the helicase activity (the hand in FigQJ could 
be decoupled (reduced) from the motor action and ATP 
intake (both remained more or less the same). In our 
terminology, these mutations involve a reduction in h 
(i.e. the overall pulling strength, Au, or the width of the 
pulse, Wh, or both), producing a reduction in the speed 
as shown in Fig. [S] We like to add that reversed motion 
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Au=1.6 



and sensitivity to the initial conditions would provide 
crucial clues on the nature of dynamics studied in this 
paper. We hope single molecular experiments in future 
would be able to probe these in detail. 



FIG. 5: Plot of c(Au,wh) as a function of the width uih of 
Eq. 10 for different values of Ait as noted against the curves. 
Solid lines are fits to c(Au, w) = cq[1 — aexp(— w/wq)]. 



of recG can also be understood in the same scheme with 
a few extra ingredients. This will be discussed elsewhere. 
Active helicases like PcrA that are involved in repair pro- 
cesses become operational when the replication process 
stalls because of, e.g., defects. Such a stalling would lead 
to a relaxation of the stalled front. The new relaxed 4>(z) 
would then act as the initial condition for the new heli- 
case recruited for repair. The sensitivity to initial condi- 
tion of a front invading an unstable phase is an important 
distinction between the metastable (pushed) and the un- 
stable cases. Whether helicases in charge of repair are 
really sensitive to and recognize these initial conditions 
need to be probed experimentally. 

To summarize, we have proposed a simple coarse 
grained model for describing helicase activity on DNA. 
The DNA is described by a local probability of zipping of 
the base pairs. The motor action of the helicase induces 
an instability around the front (Y-fork) in an otherwise 
coexisting zipped and unzipped phases. Such coexisting 
phases with an interface can in principle be created by 
pulling or by a stalled helicase on a dsDNA and therefore 
can be studied experimentally We have shown that the 
local drive created by the helicase leads to a traveling 
wave solution with a selected velocity that depends on 
the nature of the drive. Quantitative studies of forces 



[1 

[2 

[3 

[4 

[5 
[6 

ir 

[8 

[9 
[10 

[11 

[12 

[13 
[14 

[15 

[16 

[17 

[18 
[19 



* Electronic address: |email: s omenQio pb.res.in~| 
Alberts B. et al, Molecular Biology of the Cell, Garland 
Pub; 4th edition (March 2002) 

von Hippel P. H. and Delaguotte E., Cell, 104, 177-190 
(2001); Lohman T. M. and Bjornson K. P., Annu. Rev. 
Biochem. 65, 169-214 (1996). 

Patel S. S. and Picha K. M., Ann. Rev. Biochem. 69, 651 
(2000); Fass D., Bogden C. E., Berger J. M., Structure 7 
691 (1999). 

Velankar S. S. et al, Cell 97, 75 (1999); Soultanas P. etal, 
EMBO 19, 3799 (2000). 

Lohman T. M., J. Bio. Chem. 268, 2269 (1993). 
Singleton M. R. et al, Cell 107, 79 (2001). 
Bhattacharjee S. M. and Seno F., J. Phys. A 36, L181 

(2003) 

Betterton M. D., Julicher F., physics/0212045 
Chen Y. Z. etal Phys. Rev. E 56, 919 (1997). 
Bhattacharjee S. M ., J. Phys. A 33, L423(2000) 
(cond-mat/9912297); 

Marenduzzo D., Bhattacharjee S. M., Maritan A., Orlan- 
dini E. and Seno F., Phys. Rev. Lett. 88, 028102 (2002). 
Marenduzzo D., Trovato A. and Maritan A., Phys. Rev. 
E 64, 031901 (2001). 

Sebastian K. L., Phys. Rev. E62, 1128 (2000). 
Kafri Y., Mukamel D. and Peliti L., Euro. Phys. J B 
27,135 (2002). 

Lubensky D. and Nelson D. R., Phys. Rev. Lett. 85, 
1572 (2000). 

Bhattacharjee S. M. and Marenduzzo D., J. Phys. A 35, 
L349 (2002). 

Ebert U. and van Saarloos W., Physica D 146,1 (2000). 
Tripathy G. and van Saarloos W., Phys. Rev. Lett. 85, 
3556 (2000). 

Bianco P.R. etal, Nature 409, 374(2001).; Dohoney 
K.M. and Gelles J. , ibid, 370 (2001). 
Benguria R. D. and Depassier M. C. , Phys. Rev. Let. 
77, 1171 (1996). 



